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ABSTRACT: 

An  Investigation  of  elastic  buckling  of  gun-launched  motor  cases 
is  considered.  The  SATANS  finite  difference  computer  program  was 
used  for  the  analysis.  Comparative  results  from  other  buckling 
models  were  obtained  to  verify  the  results.  The  results  are  sum- 
marized in  a  set  of  design  curves  and  show  that  for  the  proposed 
design  parameters,  the  elastic  behavior  is  exceeded  before  buckling 
can  occur. 
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NOTATION 

a,  a      acceleration,  maximum  acceleration  (in  g's) 

E  Young's  modulus  of  elasticity 

G  Shear  modulus 

h  thickness  of  motor  case 

L  length  of  motor  case 

N  buckling  load  per  unit  circumference 

p.  radial  pressure  due  to  slumping  propel  1  ant 

p  escaping  breech  force  pressure 

p  motor  case  body  force  due  to  acceleration 

P  aft  edge  force 

P*  forward  edge  force 

R  radius  of  motor  case 

Wf  weight  of  structure  fore  to  motor  case 

v  Polsson's  ratio 

a  axial  stress 

a  yield  stress 

*7 

uj  density 


iii 


INTRODUCTION 

This  report  is  concerned  with  the  elastic  buckling  of  gun-launched 
motor  case  structures  which  are  subjected  to  large  axial  g  loading 
during  launch.  Acceleration  loading  up  to  10,000  g  is  typical.  The 
axial  loading  is  accompanied  by  a  lateral  loading  due  to  an  internal 
slumping  propellant  which  also  provides  additional  stiffness.  In  a 
previous  report  (ref.  1,  Appendix  D)  a  computation  of  the  buckling  load 
for  a  motor  case  with  0.08  inch  thickness,  1.5  inch  radius,  and  9.5  inch 
length  showed  buckling  occurs  at  a  stress  of  356  ksi ,  exceeding  the 
elastic  limit  of  the  high  performance  heat-treated  steel  (200-250  ksi). 
Although  that  analysis  indicated  elastic  buckling  was  not  probable,  a 
parametric  study  of  the  elastic  buckling  of  the  case  for  various  values 
of  radius,  length,  thickness  and  forward  mass  was  undertaken  in  the 
event  steels  with  higher  heat  treatment  were  used.  A  finite  difference 
computer  program,  known  as  SATANS  (ref.  2)  and  capable  of  performing 
geometrically  nonlinear  static  and  dynamic  analyses  of  arbitrarily 
loaded  shells  of  revolution,  was  used  to  compute  the  static  buckling 
loads. 


DESCRIPTION  OF  THE  MOTOR  CASE 

The  motor  case  structure  is  a  cylindrical  shell  of  4130  heat- 
treated  steel  alloy.  From  MIL-HDBK-5A,  February  8,  1966: 
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Elastic  modulus  in  tension  and  compression,  E  =  30x10  psi 
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Modulus  of  rigidity,  6  =  11x10     psi 
Poisson's  ratio  v  =  0.32 
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Density,  u>  =  0.283#/in 

The  design  parameters  furnished  by  The  Naval  Weapons  Center, 
China  Lake,  are: 

radius,  R  =  2.5  inches 

length,  L,  between  19  and  24  inches 

thickness,  h,  between  0.125  and  0.250  inches 

forward  structure  weight,  Wf,  between  45  and  70  lbs. 


DESCRIPTION  OF  THE  EXTERNAL  LOAD 


The  external  load  on  the  motor  case  structure  was  originally 
idealized  as  shown  in  the  accompanying  figure. 

Pf  is  the  force  per  unit  circum- 
ference along  the  forward  edge 
and  is  associated  with  the 
acceleration  of  the  forward 
structure. 


2*R 


m 


Po 


where  am  is  the  maximum 
m 

acceleration  in  g's. 
p    is  the  axial  distributed  force 
per  unit  circumference  due  to 
the  acceleration  of  the  motor  case  itself, 

p„     =    uiham 
rs  m 

where  oi  is  the  case  material  density. 

P.  is  the  aft  force  per  unit  circumference  obtained  by  equilibrium 
a 

considerations 

Pa    ■    Pf    +    VL 
Finally,  p.  and  p    are  the  lateral  pressures  associated  with  the 
slumping  propellant  and  escaping  breech  pressure,  respectively.     These 
two  pressures  were  neglected  in  the  analysis.     The  breech  pressure  p 
can  be  minimized  by  good  obdurator  design,  and  it  would  be  difficult 
to  determine  the  magnitude  under  any  circumstances.     The  effect  of  the 
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slumping  propel  1  ant  pressure  p^  was  examined  in  detail  and  is  the  sub- 
ject of  another  report.  It  was  determined  that  the  maximum  internal 
pressure  would  be  of  the  order  of  10  ksi  for  a  motor  case  with  a  2.5 
inch  radius,  a  10,000  g  acceleration,  and  a  20  inch  length.  Theory 
predicts  and  preliminary  analyses  verified  that  the  internal  pressure 
stiffens  the  cylinder.  This  fact,  together  with  the  uncertainty  of 
the  propel  1  ant  parameters,  provided  the  basis  for  neglecting  the 
internal  pressure  in  the  analyses. 

ANALYSIS 
The  validity  of  the  results,  which  are  given  in  the  next  section,  is 
considered  here.  The  boundary  conditions  on  both  ends  are  taken  as 
simple  supports. 

The  various  thin  shell  theories  are  associated  most  often  with 
slightly  different  approximations  and  assumptions.  For  example,  Love, 
Donnell,  Reissner,  Flugge,  Sanders  and  others  have  thin  shell  theories 
associated  with  the  particular  thin  shell  models  they  have  proposed. 
The  finite  difference  program  SATANS  used  in  this  analysis  is  based 
on  the  Sanders  model  for  thin  shells  (ref.  2),  which  includes  the  non- 
linear terms  in  the  strain-displacement  and  equilibrium  relations.  The 
results  obtained  from  SATANS  for  a  cylinder  subjected  to  an  axial  load 
are  compared  to  those  from  three  other  sources;  the  Flugge  equations 
(ref.  3),  the  classical  (linear)  buckling  equation  (ref.  4),  and 
empirical  results  (ref.  5).  The  results  for  the  critical  buckling 
load  N  and  corresponding  stress  a  are  given  in  Table  1  for  various  values 
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of  R/h.  In  all  cases,  Young's  modulus  and  Poisson's  ratio  are  30x10  psi 
and  0.32,  respectively. 


Table  1.    N  is  in  lb. /in. 


a  is  in  psi . 


SATANS 

FLUGGE 

CLASSICAL 

EMPIRICAL 

R/h  =  288 

L  =  40",  R  =  2" 

N  =  360 
a  =  52,700 

N  =  315 
a  =  45,500 

N  =  439 
a  =  64,300 

N  =  144 
a   =  20,800 

R/h  =  144 

L  =  40",  R  =  2" 

N  =  1062 
a  =  76,600 

N  =  1112 
a  =  80,400 

N  =  1755 
a  =  126,800 

N  =  633 
o  =  45,700 

R/h  =  92 

L  =  40",  R  =  2" 

N  =  3264 
a  =  149,000 

N  =  3000 
a  =  137,000 

N  =  4386 
a  =  200,000 

N  =  1656 
o  =  75,600 

R/h  =  16 

L  =  40",  R  =  2" 

N  =  46,300 
a  =  370,000 

N  =  43,300 
a  =  346,000 

Not 
Applicable 

N*  =  58,500 
a*   =  468,000 

In  all  cases  the  results  given  in  Table  1  show  the  classical  (linear) 
buckling  theory  gives  buckling  loads  which  are  not  in  agreement  with 
any  of  the  other  sources.  Its  inclusion  here  was  simply  to  demon- 
strate this  shortcoming.  It  is  also  noted  that  for  R/h  =  16  the 
classical  theory  cannot  be  used. 


The  first  three  cases  of  the  table  show  that  for  very  thin  shells,  say 
R/h  >  50,  the  SATANS  and  Flugge  theories  agree  quite  well  but  give 
buckling  loads  twice  as  large  as  the  empirical  results.  This  differ- 
ence has  been  noted  by  many  investigators  (see  ref.  5,  page  509),  and  is 
accounted  for  by  shell  imperfections  and  load  eccentricity.  For  the  motor 
case  dimensions  provided  by  NWC,  the  R/h  ratios  lie  between  10  and  20. 


Here  Flugge  and  SATANS  agree  quite  well  and  if  the  empirical  results  are 
extrapolated,  additional  agreement  is  obtained.  In  addition  to  this 
corroboration  between  theory  and  empirical  results,  it  is  noted  that  for 
these  shells  with  R/h  =  15,  the  theoretical  solution  is  on  the  safe  side 
of  the  extrapolated  empirical  result.  The  superscript  asterisks  on  the 
empirical  N  and  a   for  R/h  =  16  indicates  these  results  have  not  been 
achieved  in  the  laboratory,  rather  they  have  been  obtained  by  extrapo- 
lation to  provide  some  comparison  here. 

DESIGN  BUCKLING  CURVES 
The  results  obtained  from  SATANS  for  the  range  of  design  parameters 
given  by  NWC  are  shown  in  figures  1  through  5.  Figures  1  through  3  are 
plots  of  acceleration  versus  L/R  for  the  minimum  and  maximum  h/R  ratios 
(0.05  and  0.1  respectively).  Each  figure  is  for  a  specified  forward 
weight,  W^,  of  45,  60,  and  75  pounds  respectively.  The  buckling  curves 
are  specifically  labeled.  The  remaining  curves  denote  failure  by  yield- 
ing for  200,  250,  and  300  ksi  steels.  The  most  important  observation 
here  is  the  one  anticipated  by  the  previous  analysis;  elastic  buckling 
will  not  occur.  If  buckling  is  to  occur,  it  must  be  plastic. 

The  small  variation  of  acceleration  with  L/R  ratio  suggests  the  alternate 
plots  of  figures  4  and  5  which  show  acceleration  versus  W-  for  minimum 
and  maximum  h/R  ratios  for  fixed  L/R  ratios. 

DESIGN  TECHNIQUE 
Figures  1  through  3  or  figures  4  and  5  may  be  used  for  design  by  select- 
ing any  three  out  of  the  four  quantities,  a,  h/R,  L/R,  and  Wf,  and 
obtaining  the  remaining  fourth  quantity  from  the  lower  of  the  two  curves 
associated  with  buckling  and  elastic  yield.  In  all  cases  considered 


here  the  latter  is  the  failure  mode. 

The  values  of  h/R  equal  to  0.05  and  0.10  are  the  minimum  and  maximum 
values  of  h/R  obtained  from  the  NWC  design  parameters,  R  =  2.5  inches  and 
h  between  0.125  and  0.250  inches.     The  buckling  curve  for  h/R  equal  to 
0.075  was  determined  to  provide  a  means  of  interpolation  for  values  of 
h/R  intermediate  to  0.05  and  0.10.     The  proximity  of  the  0.075  curve  to 
the  mean  of  the  0.05  and  0.10  curves  suggests  that  linear  interpolation 
may  be  used.     As  an  example  consider  the  case  for  h/R  =  0.625,  L/R  =  7.6 
and  Wf  =  60  pounds.     From  figure  4  we  find  am  equal  to  17,600  g  and 
42,300  g  for  h/R  equal   to  0.05  and  U.lu  respectively.     3y  linear  inter- 
polation the  maximum  acceleration  for  static  buckling  when  h/R  =  0.0625, 
L/R  =  7.6,  and  Wf  =  60  pounds  is 

am    =     'Slt  I   "50       (17»600  +  42,800)  =  35,200  g. 
Because  figures  1  and  5  have  log  scale  ordinates,  direct  scaling  of  the 
curves  is  not  valid. 

The  relation  governing  the  initial  yield  curves  is 

F   ma   <"f  *  Ws>* 
0     "  A  "  A  "   2*Rh 

where  Ws  is  the  weight  of  the  motor  case,  Wf,  is  the  weight  of  the  fore 
structure,  and  the  acceleration  a  is  in  g  units.  Solving  for  the  accel- 
eration when  the  stress  is  at  yield  gives 

2*Rh  2irRh 


am  " 


wf  +  ws   ay       wf  +  2nRh  .  a.  ay 


This  equation  may  be  solved  for  any  values  of  the  design  parameters. 
It  may  be  noted  that  among  all  variables  the  maximum  acceleration, 
am,  is  linearly  proportional  to  yield  stress,  ay,  only. 
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